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We consider a two-dimensional homogeneous ensemble of cold bosonic atoms loaded inside two 
optical cavities and pumped by a far-detuned external laser field. We examine the conditions for 
these atoms to self-organize into triangular and honeycomb lattices as a result of superradiance. By 
collectively scattering the pump photons, the atoms feed the initially empty cavity modes. As a 
result, the superposition of the pump and cavity fields creates a space-periodic light-shift external 
potential and atoms self-organize into the potential wells of this optical lattice. Depending on the 
phase of the cavity fields with respect to the pump laser, these minima can either form a triangular 
or a hexagonal lattice. By numerically solving the dynamical equations of the coupled atom-cavity 
system, we have shown that the two stable atomic structures at long times are the triangular lattice 
and the honeycomb lattice with equally-populated sites. We have also studied how to drive atoms 
from one lattice structure to another by dynamically changing the phase of the cavity fields with 
respect to the pump laser. 

PACS numbers: 3T.30.-hi, 37.10.Vz, 37.10.Jk, 42.50.Pq 


I. INTRODUCTION 

Cavity QED (CQED) investigates the interaction of 
atoms with confined electromagnetic field modes. When 
atoms are coupled to a high-finesse optical resonator, the 
usual free-space dipole force they experience is strongly 
enhanced and, at the same time, the back-action of atoms 
on the confined light field cannot be ignored any longer. 
As a consequence both atoms and light must be treated 
on the same footing and the dynamics for the atomic 
motion and the cavity field becomes strongly nonlin¬ 
ear. These hybrid systems open the way to new phys¬ 
ical situations, in particular when cold atoms, bosonic or 
fermionic, are loaded inside optical cavities [lj[2]. Fur¬ 
thermore, it is possible to probe the properties of these 
systems in a non-destructive way by using the field leak¬ 
ing outside the cavity [3j. 

Trapping Bose-Einstein condensates (BEC) in laser- 
driven high-finesse optical cavities [4] has been realized 
experimentally recently [SHE]. The strong atom-cavity 
coupling enhances nonlinear effects and bistable behav¬ 
ior inn] and even chaos m set in. In the dispersive 
regime where the pump and cavity fields are both far- 
detuned from the atomic transition, the light fields im¬ 
part forces on the atoms which thus move. In turn, the 
light fields pick up phase shifts induced by the refractive 
index of this moving atomic dielectric medium. This al¬ 
ters the light forces, thus the atomic motion, thus the 
accumulated dispersive phase shifts, and this combined 
atom-field process loops self-consistently. As a result, 
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when the pump field strength is larger than some crit¬ 
ical value, the atomic cloud scatters constructively the 
pump photons into the cavity modes. This causes an 
abrupt increase of the number of photons inside the cav¬ 
ity, a phenomenon referred to as superradiance, and the 
atoms achieve self-organization into the effective opti¬ 
cal lattice created by the coherent superposition of the 
pump and cavity fields [fl4HT8] . Self-organization breaks 
the initial translation symmetry. In other words, pho¬ 
ton scattering couples the initial zero-momentum state 
of the atomic cloud to a superposition of higher recoil 
momentum states m • This effect has been used to sim¬ 
ulate the Dicke superradiance quantum phase transition 
[20H22] in a BEC-cavity system where two collective mo¬ 
tional modes of the condensate play the role of the two 
hyperfine spin states of the original Dicke model |232- 
[25]. The non-equilibrium dynamics of such systems have 
been studied in [19j 26H3Q] . Recent theoretical proposals 
have generalized the model even further by introducing 
cavity-assisted Raman coupling in order to reach phases 
such as self-organized magnetic lattices of bosons [31. or 
topologically non-trivial phases of fermions |32j [33] . 

In this article we consider a two-dimensional system 
similar to the experiment in [24] but with two crossing 
cavities in order to examine the possible formation of 
self-organized triangular and honeycomb lattices of cold 
bosons as a result of superradiance. Lately, the honey¬ 
comb lattice has attracted a lot of attention in the cold 
atom community because of its unique band structure 
mimicking Weyl-Dirac quasi-particles at the so-called 
Dirac points [34H37] . By numerically simulating the real¬ 
time dynamics of the cavity-atoms system we show that, 
depending on the relative phase (j) of the cavity fields 
with respect to the pump field, three distinct atomic lat- 
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tice structures are possible when non-interacting atoms 
and identical cavities are considered. The first one is a 
triangular lattice, the cavity fields being in phase with 
the pump laser (0 = 0). The second one is a honeycomb 
lattice with density-balanced sites, the cavities and pump 
fields being in quadrature (0 = 7r/2). The last one is a 
honeycomb lattice with density-imbalanced sites, the rel¬ 
ative phase 0 being in between zero and 7r/2. We address 
the long-time stability of these atomic lattices as well as 
the possibility of driving the atoms from one structure to 
another. 

The rest of this paper is organized as follows. In 
Sec. [II] we introduce the experimental setup that we con¬ 
sider and giv e the relevant Hamiltonian in the dispersive 
regime (Sec. |II A ). The dynamical equations of the cou¬ 
pled atoms-cavities system are derived in Sec. |IIB| They 
give rise to an effective lattice potential for the center of 
mass motion of the atoms which is studied in details in 
Sec. |II C[|IID| and |IIE| Using the symmetry properties of 
the effective lattice potential, we derive the atomic equa¬ 
tions of motion in reciprocal space in Sec. |III| The con¬ 
dition for the normal to superradiance phase transition 
is obtained by studying the linear response of the system 
in Sec. |IV| Finally, in Sec. [V| we present our numeri¬ 
cal results about the atoms and cavity fields dynamics, 
the superradiance phase transition (Sec. VA), the long¬ 
time stability of the different atomic lattices which are 
obtained (Sec. VB) as well as the possibility of switch¬ 
ing t he sy stem between these different latti ce s tructures 
(Sec. |V C) . We summarize the work in Sec. VI 


II. PHYSICAL SITUATION AND MODEL 
A. Hamiltonian 

To study the self-organization of bosonic atoms into a 
triangular or a honeycomb lattice as a result of superra¬ 
diance, we consider two high-finesse optical cavities with 
frequencies ujj (j = 1, 2) located in the xy -plane and ori¬ 
ented symmetrically about the y -axis by an angle 6 = tt/3 
(see Fig. [l]) . A two-dimensional (2D) dilute cloud of non¬ 
interacting two-level ultracold bosonic atoms (resonance 
frequency cj a , excited-state lifetime T) is then loaded in¬ 
side these cavities, which we assume do not contain any 
photons initially. A monochromatic standing-wave at fre¬ 
quency uo p (obtained by retro-reflecting a classical laser 
field propagating along unit vector y) is used to pump 
the atoms. In such a system, the atoms scatter the pump 
field photons into the cavity modes and self-organize in 
the potential wells of the effective potential created by 
the standing-wave and cavity fields. 

In the following, we assume that spontaneous emission 
processes are fully negligible so that the atomic dynam¬ 
ics is Hamiltonian. This is achieved when both the pump 
and cavity fields are far detuned from the atomic reso¬ 
nance frequency, \uj p — uo a \ T and \ujj — uj a \ T. We 
further assume that the rotating-wave approximation is 
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FIG. 1: (Color online) Schematic drawing of the setup that 
can be used to create 2D triangular and honeycomb lattices 
of bosonic cold atoms. A cloud of atoms (blue disk) is loaded 
inside two initially empty optical cavities symmetrically ori¬ 
ented by the angle 0 — 7t/3 about axis Oy. The atoms are 
then pumped by a retro-reflected classical laser field {Oy di¬ 
rection) and release photons into the modes of the cavities. 
In return, the coherent superposition of the pump and cavity 
fields creates an effective periodic potential and atoms self- 
organize themselves into the corresponding potential wells. 
A triangular or honeycomb lattice structure of atoms thus 
emerges as a result of superradiance processes. 


valid such that fast variables can be eliminated. This 
is the case when \uj p — uj a \ <C uj p and | uj p — Uj\ <C u p . 
Finally, we assume that one can adiabatically eliminate 
the atomic excited state amplitude so that atoms, ini¬ 
tially prepared in their ground state, mostly evolve in 
their ground state. This is the case when the Rabi os¬ 
cillation has a small amplitude, which happens when the 
pump detuning A a = uj p — u: a is much larger in magni¬ 
tude than the kinetic energy of the atoms in their excited 
state, the pump and cavity Rabi frequencies, and the 
pump-cavity detuning A j = uj p — ujj. Under all these as¬ 
sumptions, and since the pump is described by a classical 
field, the Hamiltonian describing the coupled dynamics 
between the cavity fields and the 2D motion of atoms, in 
their internal ground state and at center-of-mass position 
r = rx + yy , reads: 


H = 


n = 


+ 


J dr - ]T hAj a ) aj 
_A(A + A ) + ^ fV) 

2m W dy 2> A a p{ 1 
j= 1,2 

Y ^~ RF p( fi ) F j( p )( a j + a i) 

j= 1,2 a 

^X^ L Fi(f)F 2 {r){a\a 2 + a\a i), 


( 1 ) 


( 2 ) 
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FIG. 2 : (Color online) Left panel: Triangular Bravais lattice 
associated to the effective potential V e s (p, t), Eq. Q, and 
its Bravais vectors a\ and a 2 , Eq. Right panel: The 

corresponding triangular reciprocal lattice spanned by vec¬ 
tors bi and £>2, Eq. . The initial zero-momentum state of 
the atomic cloud (black disk at the center) is directly cou¬ 
pled to eighteen points of the reciprocal lattice (other colorful 
disks) by Schrodinger’s equation Eq. § . Disks with the same 
color refer to a given coupling term in Schrodinger’s equation 
Eq. Pump-cavity coupling: 771 (purple) and 772 (blue) 

terms. Inter-cavity coupling: U12 term (red). Intra-cavity 
coupling: U\ (green) and U2 (brown) terms. Pump coupling: 
Up term (orange). 


where Q p is the pump Rabi frequency, gj the atom- 
cavity coupling strength, and 4/ and ctj the atomic and 
cavity bosonic annihilation operators. The atomic op¬ 
erator is normalized to the total number of atoms AT, 
f dr 4^(r, £)T(r, t) = N. With a suitable choice of the 
origin of coordinates, the pump and cavity mode func¬ 
tions can be written as 

F p (r) = cos (k p • r) 

Fj(r) = cos (kj - r + 4>/ 2), (3) 

where 0 is a controllable phase that can be changed, for 
example, by moving the pump mirror along Oy. The 
pump and cavity wave vectors are k p = k p y and kj = 
kjbj respectively, with the unit-length vectors 

- Vs 1 

bi = sin#x + cos#y = -^-x + -y 

Vs 1 

b 2 = - sin#x + cos#y =-—x +-y. (4) 

The pump and cavity wavelengths are X p = 27r/k p = 
2ttc/uj p and A j = 27r/kj = 2 7rc/ujj, c being the speed of 
light. 


B. Dimensionless mean-field equations 

In the rest of this paper, we consider the mean-field 
regime where the field operators are replaced by their 
mean values, \h(r, t) —» (\h(r, t)) = 0(r, £) and aj(t) 

(dj(t )) = Qj(t), and quantum fluctuations are discarded. 
We next define the pump recoil energy Er = h 2 k 2 /(2m). 


Using k^ 1 , = H/Er and Er as space, time and en¬ 

ergy units, the dimensionless Schrodinger equation for 
the atomic wavefunction reads: 


id T >p(p, t ) = 


-i^ + ^ + UpCosiy 


+ Yjj= 1,2 Vj ( a j + a j ) cos V cos (bj ■ p + 4>/2) 
+Ui2(a\a.2 + aia. 2 ) cos(6^ • p + <f>/ 2) cos(&2 • p + <j>/ 2) 

+ Uj\aj\ 2 cos 2 (6'- • p + 0/2) tp(p,r) (5) 


where r = u R t, p = k p r = rx + yy, ip(p,r) = i/j(r,t)/k p , 
and b'j = kj/k p = (A p /A j)bj. Note that the normal¬ 
ization condition for the reduced atomic wavefunction is 
unchanged, J dp\(p(p,r)\ 2 = N. The various dimension¬ 
less coupling constants appearing in Eq. © are: 


n 2 p 

u ’ = ^± u ’ = 


URA a 


Vj = 


gjfl p 

-a - Ul2 

u R A a 


9i92 

URA a 


(6) 


Introducing the dimensionless cavity decay constants Kj , 
the equations of motion for the cavity field amplitudes 
read: 


id r ai = —(Si + iK\)ai + NUi 2 h 2 a 2 + Nrjil lp 
id r a 2 = ~(S 2 + iK 2 )a 2 + NUi 2 h 2 ai + Nr] 2 I 2p , (7) 

where Sj = A j — NUjlj is the shifted cavity resonance 
frequency and where: 


h = 

h 

[ dp \ip(p,r)\ 2 cos 2 (fy -p + 0/2) 


h2 = 

h 

[ dp \ip(p,r)\ 2 cos(b\ -p + ct>/2) cos(b' 2 

■ p + </>/2) 

IjP = 


1 dp \<p(p,T)\ 2 cosycosibj-p + <j>/2). 

( 8 ) 


Note that, because of the normalization of the atomic 
wave function, the /-integrals do not depend on the ac¬ 
tual number of atoms N. 


C. Effective potential for atoms 

From Eq. (J5| , one immediately sees that the pump and 
cavity fields create an effective potential for the atom 
center-of-mass dynamics: 


Ves(p,T) = Up COS 2 jj + J2j=l,2 Vj\ a j | 2 COS 2 ( 6 ' ■ p + fi/2) 

-\-U \ 2 (ol\ol 2 + ol^oli) cos^i • p + 0/2) cos(b' 2 • p + 0/2) 

+ Ej=i ,2 7 7j( Q q + aj)cosycos(b'j • p + 0/2). (9) 

This effective potential can be recast under the simpler 
and suggestive form 


V e s(p, r) 


h\Q(r, t)\ 2 

ErAcl 


( 10 ) 
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featuring the complex Rabi field amplitude 

D(r,t) = Q p Fp(r) + E (ii) 

i=*i»2 

associated, in the mean-field regime, to the coherent su¬ 
perposition of the pump and cavity classical standing- 
wave electrical field amplitudes. As one can see, the 
pump and cavity fields are balanced when rij = |ay| 2 = 
= U P /U j {j = 1,2). 

We will restrict our subsequent analysis to a red- 
detuned pump field, A a < 0, for which the minima of 
the effective potential, Eqs. ([9 ])-([Tq]) , are found at the 
maxima of the norm of the Rabi field, Eq. ( pT| ), and will 
thus be the deepest. 


D. Lattice symmetry of the effective potential 


Let us define the two vectors a'- = (Xj / X p )aj (j = 1, 2) 
with 


“i = »<—»+ 


sinO ' 


= + < 12 > 


and \dj | = 47r/\/3. As readily checked, they satisfy the 
property a[ • b'- = 2^5^. As a consequence, and by 
construction, the cavity mode functions Fj are left in¬ 
variant when p is translated by any integer combination 
of the o!y The pump mode function F p is also left in¬ 
variant provided a'-.y is a multiple integer of 2i r. The 
simplest choice, which also fulfills the conditions lead¬ 
ing to Eq. ©. is Xj = X p . If this condition is not strictly 
obeyed, F p will still be approximately invariant, to a very 
good accuracy, if \k p — kj\L 1, where L is the size of 
the atomic cloud. In the rest of the paper we will as¬ 
sume this condition to hold such that the approximation 
Ai ~ A 2 ~ Xp = X is perfectly justified. In turn, it means 
that the cavity and pump fields have almost the same 
frequencies, \Aj\ <C uj p [38]. 

Within these assumptions, it follows that the Bravais 
lattice associated to the effective potential Wff(p, r) is 
C = {R = J2j see Fig. [2J It is triangular 

and its unit cell is A = {JT V A] 0 < Vi < 1} with dimen¬ 
sionless area S = 2(27r) 2 /\/3. The corresponding recipro¬ 
cal lattice 7Z = {K = X] j Pj bj 5 Pj €= is also triangular 
and its first Brillouin zone is 13 = \ u j I ^ 1/2} 

with dimensionless area E = (27 t) 2 /S = y/3/2. 


E. Case of identical cavities 

In the following, we will examine the simplest case 
where the cavities have exactly the same characteristics, 
lji = UJ2 = cj c , gi = p2 = Qc an d ki = ^2 = ft c - Then a 


reflection symmetry about axis Oy (x —>> —x) amounts to 
exchange the cavity fields, Aq Aq. If we further assume 
that the cavities are initially empty and the initial atomic 
cloud is also symmetric with respect to Oy , then the sub¬ 
sequent time evolution will always enforce aq = a 2 at all 
times [39,. In this case the Rabi field, Eq. ( [TT] ) , writes: 


f }(p, r) = Vtp cos y + 2f 1 c (r) cos( 


y + (t> 
2 



) (13) 


with Q c (t) = g c OLi{r) = < 7^2 (t). It is immediately 
seen that, in this case, the effective potential is al¬ 
ways reflection-symmetric about planes located at x = 
277 / 77/^/3 (n G Z) and in particular about axis Oy (x —>> 
—x). It is also easy to prove that the potential for 0 + nir 
(' n G Z) is obtained by a mere translation of the po¬ 
tential for (j). Furthermore, the potential for —0 is sim¬ 
ply obtained by reflecting the potential for 0 about axis 
Ox (y —» —y). This means that, for all practical pur¬ 
poses, the range of cavity phases can be restricted to 
0 < 0 < tt/2. 

For (j) = 7t/2, the Bravais unit cell hosts two energy- 
balanced minima which can be labeled A and 5, so that 
the full lattice of minima is now a honeycomb lattice 
made of two shifted A and B triangular lattices, see 
Fig- i Interestingly enough, this honeycomb structure 
does not depend on the relative weight between the pump 
and cavity Rabi amplitudes [40] . Indeed, when 0 = 7t/2, 
the potential is symmetric under p —)> 7 ry — p. This im¬ 
mediately ensures the existence of an even number of 
energy-balanced minima in the Bravais unit cell, here 
two. 

When 0 is continuously decreased from 7t/2, the A 
and B minima become energy-imbalanced, still retaining 
a nice-looking honeycomb structure when cj> is not too 
small. This situation is interesting for producing bands 
with non-vanishing Chern numbers. When <p gets closer 
to zero, the honeycomb structure is lost for all practical 
purposes and the lattice of deepest minima is triangular. 

Irrespective of the value of 0, and for a priori differ¬ 
ent mean-field values ay, it is worth mentioning that the 
effective potential gets shifted by d\/2 when the sign of 
0 L 2 is flipped, by < 22/2 when the sign of aq is flipped and 
by (ai + ^ 2 )/2 = 27ry (which corresponds to a shift of A 
along axis Oy) when both signs are flipped. It is straight¬ 
forward to check that Eqs © and 0 indeed remain in¬ 
variant when the sign of aq and/or a 2 is flipped and the 
corresponding translation by aj/2 is implemented on the 
mode functions and atomic wave function (/?(p, r). There¬ 
fore, for each set of coupling parameters, there are four 
possible solutions associated to given numbers of cavity 
photons {rij = |ay| 2 ), all related by translations along 
ai/2 and/or 0 , 2 / 2 . Depending on the initial conditions, 
the system may select any one of these solutions. 
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FIG. 3: (Color online) Effective potential V e a (p, r) given by Eq. § and its lattice structure of minima (in units of the common 
wavelength A associated to the cavity and pump fields) for three values of the cavity phase 4> appearing in the mode functions 
Fj(f) (see text). Right panels: for (j) = 7r/2 a honeycomb structure is obtained for the potential minima which have all same 
depths. Middle panels: for lower values, e. g. 0 = 7t/4, a honeycomb structure is obtained but with energy-imbalanced minima. 
Left panels: for <f) = 0, a triangular lattice is obtained for the potential minima. The dimensionless coupling parameters are 
U p = -4, Uj = U 12 = -0.02, r]j = -y/2/b, see Eqs ©> and dj (r) = (a*(r)) = 15 (j = 1,2). With these values, the pump and 
cavity fields are almost balanced since U p /Uj = 200 and rij — \otjf : = 225, see text after Eq. ([ID- 


iii. SCHRODINGER EQUATION IN 
RECIPROCAL SPACE 

We now rewrite Schrodinger’s equation, Eq. (§, by 
expanding the atomic wave function <p(p, r) in reciprocal 
space: 


It is straightforward to see that q is actually con¬ 
served during the evolution. Indeed, the effective poten¬ 
tial T4ff(p, r) being periodic under Bravais translations, 
it cannot scatter and change q. Substituting Eq. (14) 
into Eq. © and using Eq. we obtain the following 
dynamical equations for the coefficients Cg(q,r): 


<p(p,t) = y E E C R (q,T)e l ^+ K >P (14) 

where N [4 is the number of atoms in the unit Bravais cell 
A. The normalization conditions now read: 

f dp\ip(p,T)\ 2 = N a (15) 

Ja 

E E KW,t)| 2 = 1. (16) 

K^TiqeB 


U n 


+ C K- 2b 1 -2b 2 ) + E 1 4 ~ ( 2C K + e#C 'x_2b i + e l ‘ t ’ C K+2b j ) 


id T Cg — (q + K) 2 Cg + J + Cg +2 y )i+2b2 ~r 2b 1 -2b 2/ 

_l_ a j ) ( i<f >/^ A A 1 . 1 P i( t>/ 2 n ^ . 1 ^ , .A 

^ 4 V ^K-bi-bs-b, ^ e ^K+bi+bs+bj ^ e ^K+bi+bs-b, ^ e ^K-bi-ba+bj-y 


3 = 1,2 
+ < 


3 = 1,2 


t/i 2 (a*a 2 +0:2^1) / +C f K+b b) 

4 y K-bi-b 2 K +bi+b 2 K-bi+b 2 rC+bi —b 2 y 


( 17 ) 
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where the g-dependence of the C-coefficients is the same 
on both sides of the equation. The dynamical equations 


J 


for the cavity fields aq, Eq. ([7]), remain the same but with 
the integrals in Eq. wh now reading: 


1 = - 
J 4 


E ^( 2 Cg + e-VCg^+eVC 


K- 2b 


) 


Kq 


(18) 


t+b 2 ^ e bi-b 2 


^K+ bi-b 2 + ^K-bi+bs) 


^ bH 

K4 

I 0 P = 4 ^^K+bi+bs+bj- + e t4>/2 c^_ bi_b 2 _b j + e * </>/2 C f ^_^ i _^ 2+ ^. + ^^C'x+bi+bs-b,) 




r 


Generally, in superradiance and self-organization prob¬ 
lems, the atomic system is initially homogeneous (zero- 
momentum state), meaning that only the C-coefhcient 
with Km 0 and qm 0 is nonzero initially (and equal to 
1 because of the normalization condition). In the course 
of time, as a result of superradiance, higher-momentum 
states will be occupied. For our system, the dynamical 
equations ( |17| ) couple Cg =Q (q = 0, t m 0) to the eighteen 
points of the reciprocal lattice 1Z shown in Fig. [2] Each 
of these points will be in turn coupled to other points 
in the reciprocal lattice with even higher momenta. In 
the following sections, we restrict our calculations to mo¬ 
menta transfers less than 3 hk, i.e. to reciprocal lattice 
points K = ^ZjPjbj with \pj\ < 3 and \pi — P 2 I < 3. 
These thirty-seven points in total include the initial zero- 
momentum state, the eighteen states to which this state 
is directly coupled, and the eighteen other states to which 
the states with momentum hk (the points closest to the 
initial zero-momentum state) are directly coupled. 


to destabilize the initial state for given cavity parame¬ 
ters. For this purpose we consider the linear response 
of the system to perturbations added to the initial cav¬ 
ity fields and atomic wave function. We thus substitute 


a 


(o) 


a 


(o) 


+ 5 aj and C 


(o) 


j ' -J —J - -_K ■ — K into the d y- 

namical equations Eqs Q and (17) and in the integrals 


C l ?+SCg 


given in Eq. (18). Our initial state is a^ = 0 (empty 
cavities) and C^ = 5^5^ (homogeneous cloud). For 
this initial state, Ij = 1/2 and 1 12 = 0. Since q cannot 
change under evolution, as explained in Sec.|III[ its value 
remains zero and will thus be omitted. The linearized 
perturbed equations then read: 


IV. INSTABILITY AND PHASE TRANSITION 
CONDITION 

In this section, we study the conditions for superradi- id r 5aj 

ance to take place and focus on the pump strength needed with 


= — A,— 


NUn 


■ IK 


i'jSotj + NpjSIjp, (19) 


Slip = ~ 


+ 


4 ETk ( e l ^ /2 S C K+b 1 +b 2 +b j Jrel<t>/2 SC i l_y }i _y }2 _y )j + e ^SCg^ _£ 2 ^ + C^SC + £ 2 ) 

K 

1 


- V SC# (e-^C^ - - - + e^ /2 C ( - 0) ~ - - + e~ i<l>/2 C ( ? - - - + e i<t>/2 C ( ?~ ■ . Y (20) 

4 ' K'+bi+b2+bj K — bi—b2—bj K — bi—b2 _ l _ bj K- |-bi-|-b2 — b jj 

K 

tnd 

= ( g2 + Y )ic t + T( ic M^ + ic *-^-^) + E 

3 = 1,2 

+ e~ i(f)/ 2 C { ? . . . + e i(t)/ 2 C { ? . - . + e -^/ 2 cl 0) . . - ) ( 21 ) 

K+bi+bs+bj K+bi+bs-b,- K-b 1 -b 2 +b j J V > 
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FIG. 4: Imaginary (top panels) and real parts (bottom panels) of the eigenvalues uj of stability matrix M as functions of 
the (negative) dimensionless pump coupling parameter U p (red-detuned pump). We have used thirty-seven reciprocal lattice 
vectors to compute and diagonalize M, see discussion end of Sec. |III| Positive imaginary parts signal instability of our initial 
state (empty cavities, homogeneous atomic cloud). The dimensionless cavity parameters are Uj = U12 = —1.74 x 10 -2 , 
A j = —4.7 x 10 3 and Kj — 480 (j — 1,2). The total number of atoms is N = 1.5 x 10 5 . 


Eqs (20) and (21) only involve reciprocal lattice vectors 


of the form K + Q with a transfer lattice vector Q = 
Qjbj satisfying Qj = 0 , ± 1,±2 and (Q 1 - Q 2 ) = 
0, ±1. Then it can be seen that all phase factors can be 
gauged away: 


g 


( 0 ) 

K+Q 


^c' 


( 0 ) 

K+Q 


— e ir(Q)<f >/2 q(0) 


K+Q 


SC- _V sc- - — e* r (<W / 2 hC- - 

°^K+Q ^ °^K+Q e 0 Ly K+Q’ 


( 22 ) 

(23) 


where (Q 1 +Q 2 ) = 4n + r(Q), with n G Z, |r(Q)| < 3 and 
r(Q) the (positive or negative) remainder with smallest 
absolute value. In other words, the stability of our initial 
state is independent from the cavity phase which can thus 
be set to <p = 0 in the stability analysis. 

Let us define the column vector Y = ({SCg}, {£aq}) T , 
where the superscript T denotes transposition. The lin¬ 
earized equations Eqs ( p~9| ) and ( | 2 l| ) can then be rewritten 
in compact form id T X_ = MX where the stability matrix 
M controls the dynamics of the (doubled) perturbation 
column vector X = (Y, Y*) T around our initial condi¬ 
tion. As one can see, the perturbation will grow, and the 
initial system is unstable, as soon as the eigenvalues uj of 
M have a positive imaginary part. The matrix M has 
the following 4-block structure 


M = 


F 

-G* 


G 

—F* 


(24) 


where the star stands for complex conjugation and where 
the sub-matrix G is symmetric G T = G. It is easy to 
show that cr x M(Jx = — M* , where the usual entries of 
the Pauli matrix a x have been replaced by the null and 
identity matrices. As a consequence, uj and — uj* are both 


eigenvalues of M, with eigenstates X LC and cr x A*, and 
the real parts must come in opposite pairs. 

Fig. [4] shows the real and imaginary parts of the eigen¬ 
values of the stability matrix M as functions of the pump 
coupling parameter U p and obtained for fixed typical cav¬ 
ity experimental parameters. Note that we have not re¬ 
produced in these plots eigenvalues with large imaginary 
parts (equal to the cavity damping rates Kj). As ex¬ 
pected from the symmetry of Af, one can see that the 
real parts come indeed in opposite pairs. However, con¬ 
trary to what a quick and misleading glance at the plots 
may let think, imaginary parts should not come in oppo¬ 
site pairs. Indeed, the sub-matrix F becomes Hermitian, 
= F, only when the cavities are lossless which implies 
the additional symmetry a z A4a z = . In this partic¬ 

ular case, both uj and uj * are eigenvalues and the imagi¬ 
nary parts would also come in opposite pairs. This is not 
the case here since the cavities are lossy (kj 7 ^ 0 ). And 
indeed, a careful check shows that seemingly opposite 
values of the imaginary parts are in fact slightly differ¬ 
ent. We have checked that this difference approaches zero 
when the cavity damping rates Kj go to zero and the cavi¬ 
ties become lossless. Another important point to mention 
is that all the imaginary parts (positive or negative) are 
either two- or four-fold degenerate. The two-fold degen¬ 
eracy is always present and comes from the fact that uj 
and — uj* are both eigenvalues. The occasional four-fold 
degeneracy comes from an additional spatial symmetry 
of the system: when cavities are identical, the system 
is reflection-symmetric about axis Oy regardless of the 
value of 0, see discussion after Eq. (13). Therefore, in 
some cases, there are two additional eigenstates of M co¬ 
existing with X LC and and sharing the same imag¬ 

inary part. They are linear combination of the images of 





















these eigenstates obtained by reflecting them about axis 
Oy. Writing K = ^ZjPjbj, these images are obtained by 
exchanging pi P 2 and Sa± Sa 2 . We have confirmed 
this fact by verifying that the 4-dimensional eigenspace 
associated to four-fold degeneracies is indeed stable un¬ 
der reflection about axis Oy. 

For the cavity parameters used in the plots, the top- 
left panel of Fig. [4] suggests that our initial state is stable 
and superradiance cannot take place when \U P \ < 1. A 
few other stable regions appear when \U P \ increases fur¬ 
ther. For larger values \U P \ > 15, a continuous instability 
develops, see top-right panel of Fig. [4] where eigenvalues 
with growing positive imaginary parts are visible. The 
threshold for unstable behavior at weak pump fields does 
not change when the number of reciprocal lattice vec¬ 
tors used to compute and diagonalize M. is increased. 
However, by expanding the momentum transfer limits 
and including more reciprocal lattice vectors, additional 
eigenvalues with positive imaginary parts appear in the 
range \U P \ > 1. 

A word of caution about the conclusions drawn from 
the linear stability analysis is necessary. Indeed, the pic¬ 
ture may change when subsequent nonlinear terms are 
included to refine the analysis. As a matter of fact, as 
we have numerically checked, our initial state is still dy¬ 
namically unstable when \U P \ is chosen in between the 
linear instability regions visible in the top-left panel of 
Fig. i However, when \U P \ is chosen within the linear 
instability regions, a better atomic lattice structure is 
obtained numerically. 


V. SELF-ORGANIZED TRIANGULAR AND 
HONEYCOMB LATTICES 


We present now the results of our numerical simula¬ 
tions of the dynamical equations ([7| and © with the 


integrals given in Eq. (18). As mentioned in Sec. Ill 


we only consider momenta transfer to atoms up to 3 hk, 
which amounts to include thirty-seven points of the re¬ 
ciprocal lattice in our simulations. Starting from empty 
cavities and a homogeneous atomic cloud, the system un¬ 
dergoes, for appropriate laser fields strengths, a super¬ 
radiant phase transition and subsequently self-organizes 
into either a triangular or a honeycomb lattice of atoms 
depending on the value of the cavity phase 0. 


a sharp and stable lattice of atoms. This is because the 
pump field is too weak and cannot feed a sufficient num¬ 
ber of photons inside the cavities. This problem is over¬ 
come by decreasing further the pump field strength down 
to U p = —3 see Fig. [5] As shown in the top-left panel 
and the insets, the number of photons inside the cavities 
increases steadily and reaches rij ~ 10 after a short time 
r^8. The atomic density has been computed at three 
different times. At small times, and before superradi¬ 
ance takes place (r = 2), the atomic cloud self-organizes 
according to the minima of the pump standing-wave po¬ 
tential alone. As superradiance takes place and the num¬ 
ber of photons inside the cavities increases (r = 25 and 
r = 130), the superposition of the pump and cavity fields 
creates a sufficiently deep potential and the atoms self- 
organize into the triangular lattice of potential minima, 
see left panel of Fig. [ 3 ] It should be noted that, within the 
time span shown in Fig. [5j the number of cavity photons 
fluctuates a lot, even if it remains sizable, and the effec¬ 
tive potential fluctuates too. This is because the system 
has not yet reached the steady state. However these fluc¬ 
tuations do not alter the triangular nature of the lattice 
of potential minima and, in turn, the triangular nature 
of the self-organized atomic lattice. 


2. Cavity phase 0 = 7r/4 


As suggested by the results of the linear stability anal¬ 
ysis developed in Sec. |IV[ superradiance takes place ir¬ 
respective of the value of the cavity phase 0. However 
the lattice structure of the potential minima is triangu¬ 
lar only when 0 = 0, otherwise it is rather hexagonal 
(see Sec. |Iie| . To check this point, we have considered 
the same initial state and we have numerically solved 
Eqs. Q, @ and ( fl8| ) with 0 = 7r/4. In this case, we 
expect superradiance to organize atoms according to the 
energy-imbalanced wells of the honeycomb potential, see 
middle panel of Fig.[3]for an example of such a potential, 
with more atoms trapped in the deeper wells. This is in¬ 
deed the behavior observed in Fig. [6] where the obtained 
atomic density distributions are shown at times r = 25 
and r = 130 and display the form of a density-imbalanced 
honeycomb lattice. 


3. Cavity phase 0 = tv/2 


A. Transition to the superradiant state 

1 . Cavity phase 0 = 0 

To check that the system does enter a superradiant 
phase, we monitor the dynamics of the system for a pump 
field strength larger than the lower bound found in the 
linear stability analysis, see Sec.|IV| Our numerical simu¬ 
lations confirm the superradiant state for U p ~ —1. How¬ 
ever the effective potential is not deep enough to support 


By the same token, when 0 = 7r/2, the wells of the 
honeycomb potential have all the same depth and we 
expect superradiance to drive the atoms into a density- 
balanced honeycomb lattice, see right panel of Fig. [3| 
This is indeed what our numerical results show, see Fig. [7] 
where the atomic distributions obtained at times r = 25 
and r = 125 are given. 

One may have noticed that, for the same coupling pa¬ 
rameters, the atomic honeycomb lattices have lower con¬ 
trasts than the triangular one. As one can see in Fig. [3J 
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FIG. 5: (Color online) As a result of superradiance, an initial homogeneous atomic cloud loaded inside empty identical cavities 
self-organizes into a triangular lattice when the cavity phase is set to 0 = 0. The pump dimensionless coupling strength is 
Up — — 3 and all other parameters are the same as in Fig. [4] Atoms scatter pump photons into the cavity modes and, after 
a short while, superradiance and self-organization take place. Top-left panel: cavity photon numbers rij = |a^| 2 (j = 1,2) 
as functions of time r = WRt. The insets show the logarithms of m and 712 as a function of time and the horizontal grid 
lines mark the values 1, 10 and 100. From a mathematical point of view, the symmetry of the dynamical equations and of 
the initial state enforces m = n 2 at all times as confirmed here by the numerical simulation. Top-right panel: atomic density 
distribution at time r — 2, before superradiance takes place. The atomic cloud self-organizes according to the minima of the 
pump standing-wave potential alone since the cavity fields are almost zero. Bottom-left (r = 25) and bottom-right (r = 130) 
panels: atomic density distributions, and their corresponding color codes, after superradiance and self-organization have taken 
place. The atomic lattice, just like the lattice of minima of the effective potential, is clearly triangular. 



FIG. 6: (Color online) Same as Fig. [ 5 ] but with a cavity phase set to 4> — 7t/4. Left panels: cavity photon numbers rij — \aj\ 2 
(j = 1, 2) as functions of time r = ujRt. Middle (r = 25) and right (r = 130) panels: after superradiance and self-organization 
have taken place, an atomic honeycomb lattice with density-imbalanced consecutive sites is formed. 
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the honeycomb minima are already shallower than the 
triangular ones when the number of cavity photons is the 
same. Here the dynamics develops less cavity photons for 
the honeycomb lattice than for the triangular one. More¬ 
over, for the honeycomb lattice, the atoms are distributed 
over twice many sites than for the triangular lattice. All 
these reasons conspire to produce less contrasted atomic 
lattices for the honeycomb structure than for the trian¬ 
gular one. For the parameters chosen here, we infer from 
Eq. (IT) that the pump and cavity fields have same or¬ 
der of magnitude when rij = U p /Uj ~ 172. As seen from 
Figs. [5][7| this is the case for 0 = 0 where rij fluctuates 
roughly around 175. For 0 = 7 r /2 and 0 = 7 r/ 4 , rij fluc¬ 
tuates roughly around 50, or a bit less, and the pump 
field is larger than the cavity fields. 


B. Long-time stability of the different lattice 
structures 


As seen in Figs. |5][7} the cavity photon numbers rij suf¬ 
fer strong temporm fluctuations. This means that the ef¬ 
fective potential fluctuates too and thus, in turn, so does 
the atomic distribution. In fact, the system needs more 
time to reach the steady state. Fig.[ 8 ]shows the long-time 
evolution of the system when 0 = 0 (triangular lattice). 
The steady state is reached around time r ~ 2000 with 
a number of cavity photons ni = ~ 2800. In this 

case the cavity fields largely dominate over the pump 
field and the wells of the effective potential organize in a 
rectangular array rendering the (exact) triangular sym¬ 
metry less apparent as is shown in Fig. [ 8 ] However, the 
atoms accumulate in the deepest minima of this effective 
potential which still form a nice regular triangular lat¬ 
tice. The contour plots in Fig. [ 8 ] show an example of the 
atomic density distribution and the effective potential at 
time r = 3000. This distribution does not change in time 
once the system has reached the steady state. 

When 0 < 0 < 7r/2, as we have seen, atomic honey¬ 
comb lattices with density-imbalanced consecutive sites 
are synthesized. However these structures are not sta¬ 
ble in the long-time limit and atoms re-organize them¬ 
selves into the deepest potential wells which form a tri¬ 
angular lattice. As an example, Fig. [9] shows the atomic 
density distribution and the effective potential obtained 
at time r = 1500 for 0 = 7r/4. As clearly seen the 
density-imbalanced atomic honeycomb lattice found at 
early times, see Fig. [ 6 j has disappeared in favor of the 
triangular atomic sub-lattice of deepest wells. 

These density-imbalanced atomic honeycomb lattices 
destabilize in favor of the triangular lattice after a certain 
latency time T. After the time T, there is no visible occu¬ 
pation of the shallower lattice sites |4l j. The atomic lat¬ 
tice is purely triangular and becomes the stable structure 
in the long-time regime. The latency time T depends on 
the potential energy difference between consecutive sites. 
Since this potential mismatch decreases when 0 increases, 
T gets longer when 0 —>► 7 r/ 2 , see Fig. 10 Note however 


that the cavity photon numbers take a time longer than 
T to reach their steady-state. 

When 0 = 7 r/ 2 , the effective potential is hexagonal 
with perfectly energy-balanced minima. We thus expect 
a stable density-balanced atomic honeycomb lattice to 
emerge from the self-organization process. As one can 
see from Fig. ED this is indeed the case, but the cav¬ 
ity photon numbers now reach their steady-state after 
a much longer time. They even fluctuate a lot dur¬ 
ing their temporal evolution and induce in turn fluctu¬ 
ations of the atomic lattice. These atomic density fluc¬ 
tuations get strongly reduced after the time r = 6000 
but, even at time r = 10 4 , the system has not yet fully 
reached its steady-state. As one may have noticed, our 
numerical results produce different cavity photon num¬ 
bers n i 7 ^ ri 2 even if their average behavior is the same 
and their asymptotic values are equal. This asymmetry 
is due to the sensitivity of this system to numerical er¬ 
rors. However, the difference in m and n 2 approaches 
zero as the system reaches its steady state. 


C. Driving atoms from one lattice structure to 
another 

So far, we have considered the dynamical behavior of 
the cavity-atom system for a fixed pump strength, cho¬ 
sen above the threshold for superradiance, and we have 
studied the self-organization of atoms into triangular or 
hexagonal lattice structures for different cavity phases. 
We now investigate the possibility of driving atoms from 
one lattice structure to another by changing in time both 
the pump strength and the cavity phase. The protocol 
we explore is the following. The pump strength increases 
linearly in time from 0 to some final value Uf over a cer¬ 
tain time window 3At. Meanwhile, the cavity phase is 
set at 0i initially. After superradiance takes place and 
enough number of photons are produced inside the cavi¬ 
ties, the phase is then changed linearly in a time interval 
of At to its final value 0/. Finally, the cavity phase keeps 
this target value for the rest of the time sequence. 

We present below the time evolution of the cavity pho¬ 
ton numbers and the atomic distribution obtained nu¬ 
merically for this experimental protocol for 3 possible 
cavity phase values: 0 = 0, 0 = 7 r /2 and some interme¬ 
diate value that we choose to be 0 = 7 r/ 4 . 


1. Switching the cavity phase between 0 and 7r/4 


When the system starts with 0^ = 0, the atoms or¬ 
ganize in a triangular lattice. By switching the cav¬ 
ity phase to 0/ = 7 r/ 4 , a new optical potential is pro¬ 
duced featuring addi tional shallower wells. However we 
know from Sec. VB that when 0 < 0 < 7 r/ 2 , the stable 
structure at long times is the triangular lattice formed 
by atoms occupying the deepest potential wells. Since 
atoms have already self-organized in a triangular lattice, 




11 



FIG. 7: (Color online) Same as Fig. [H] but with a cavity phase set to 0 = tt/2. Left panels: cavity photon numbers rij = \aj \ 2 
(.7 = 1, 2) as functions of time r = ujRt. Middle (r = 25) and right (r =a 125) panels: after superradiance and self-organization 
have taken place, an atomic honeycomb lattice with density-balanced consecutive sites is formed. 
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FIG. 8: (Color online) Long-time dynamics of the system shown in Fig. [H] (^ = 0). Left panel: cavity photon numbers rij = ]aj | 2 
(j = 1, 2) as functions of time r = ujRt. The cavity steady-states are reached around a time r — 2000. Middle panel: contour 
plot of the (stable) atomic density distribution obtained at time r = 3000. Right panel: effective lattice potential V e & at 


r = 3000. 


changing the cavity phase from 0 to 7 t /4 will not affect 
the atomic distribution which thus remains triangular. 
On the other hand, when the system starts instead with 
c t>i = 7 t /4 and the cavity phase is changed to <fif = 0, 
the transition to the triangular lattice is inevitable and 
is illustrated in Fig. [l2j As the cavity phase decreases, 
the secondary minima of the effective potential obtained 
for rf>i = 7 t /4 become shallower and shallower and are 
occupied by less and less atoms. The triangular lattice 
becomes the prominent structure around r « 280, well 
before (j> reaches the target value 0, which triggers a con¬ 
siderable rise of the cavity photon numbers. 

2. Switching the cavity phase between 0 and tt/2 

In this case, both the triangular and density-balanced 
honeycomb lattices are stable structures of the static sit¬ 
uation. Fig. p~3] shows what happens when we start from 
the triangular lattice (</>* = 0) and gradually increase 
the cavity phase to <fif = tt/2. As one expects from the 
previous case, the atomic lattice structure remains tri¬ 
angular as long as the cavity phase has not reached its 


target value, 0 < (j)(r) < tt/2. Interestingly enough, once 
the cavity phase has reached 0/ = 7r/2, the atomic den¬ 
sity distribution starts to oscillate between the wells of 
the two triangular sub-lattices of the energy-balanced ef¬ 
fective potential, taking the form of density-imbalanced 
atomic honeycomb lattices in between. These two tri¬ 
angular sub-lattices are simply shifted by — (a\ + c?2)/3. 
As a net result, the atomic density looks like ” blinking” 
between two shifted triangular lattices, a sign of a bi¬ 
stable behavior. This blinking does not stop if the pump 
strength is kept fixed, rather than linearly increasing, 
after the cavity phase has reached <fif = tt/2. It even 
speeds up if the pump strength is ramped up faster. The 
bottom panels in Fig. [I3| display the atomic densities at 
times r = 402 and r = 404, showing that, for the ex¬ 
perimental parameters chosen, the oscillation period be¬ 
tween the two triangular lattices is comparable to the 
pump recoil time . It is worth recalling that, once 
the cavity phase is set to 7r/2, the minima of the effec¬ 
tive potential have same depth and display a honeycomb 
structure, irrespective of the cavity photon numbers (as 
long as they are not zero). This case should not thus be 
confused with the previous case (0 < 0 < tt/2) where 
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FIG. 9: (Color online) Long-time dynamics of the system shown in Fig. [6] (0 = 7 t/ 4). Left panel: cavity photon numbers 
n j ~ \ a j | 2 {j ~ B 2) as functions of time r = cjRt. The cavity steady-states are reached around a time r = 4000. Middle panel: 
contour plot of the atomic density distribution obtained at time r — 1500. Right panel: V e s at r — 1500. As one can see, the 
density-imbalanced atomic honeycomb lattice obtained at short times in Fig. [6] is not a stable structure at long times. The 
atoms re-organize in the deepest potential wells of the effective potential shown in the right panel and form a stable triangular 
lattice. 
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FIG. 10: Dimensionless latency time T against the cavity 
phase 4>. After T, atoms have re-organized into the deep¬ 
est wells of the energy-imbalanced effective potential obtained 
for 0 G ]0, 7r/2[ and exhibit a triangular lattice structure. This 
time T increases when the energy mismatch between the prin¬ 
cipal and secondary potential minima decreases, i.e. when (f> 
increases. It diverges when 0 —7r/2. The atomic honeycomb 
lattice obtained at 4> = 7r/2 is a stable structure. 


the lattice potential has minima with different depths. 
Our present results then show that, for time-driven pa¬ 
rameters, the density-balanced honeycomb structure is a 
transient state between two stable triangular structures. 

Fig. [14] shows what trivially happens for the reverse 
process, i. e. when we start with = it/2 and decrease 
gradually the cavity phase to <fif = 0. Once the cav¬ 
ity phase departs from 7t/2, the effective potential has 
minima with different depths. As long as the potential 
mismatch is weak enough, the initial density-balanced 
honeycomb structure resists but, with increasing poten¬ 
tial mismatch, finally gives in and a density-imbalanced 
honeycomb lattice is formed around r = 250. Then, as 
the cavity phase decreases further, the depth of the sec¬ 
ondary minima, as well as the number of atoms sitting on 
them, decreases rapidly. Before the cavity phase finally 
reaches 0 = 0 at r = 300, the atoms choose the deeper 
triangular sub-lattice which becomes the stable structure 


as expected. 


3. Switching the cavity phase between 7 t /4 and 7 r /2 


We address now the last case, namely the transition 
between the two possible, density-balanced or density- 
imbalanced, atomic honeycomb lattices. We first con¬ 
sider switching the cavity phase from fa = 7r/4 to 
c/)f = 7t/ 2. In this case, some of the atoms are sitting in 
the shallower potential wells of the effective honeycomb 
potential while the majority of them are hosted in the 
deeper ones. By increasing the cavity phase, the shal¬ 
lower sites become deeper. However, the atoms in the 
deeper wells tend to stay where they are (see Sec. V C 1) 
until the cavity phase becomes exactly tt/2 and all po¬ 
tential wells have same depth. At this point, the atomic 
lattice is hexagonal with almost density-balanced sites. 
Then, for the rest of the time sequence, atoms start os¬ 
cillating in time between the two stable triangular sub¬ 
lattices, featuring the ’’lattice blinking” already observed 
in the preceding Sec. |V C 2| 

When the system starts with fa = 7r/2, the initial 
atomic lattice is a density-balanced honeycomb lattice. 
By gradually decreasing the cavity phase, the poten¬ 
tial wells of one of the two triangular sub-lattices be¬ 
come shallower. However, the initial density-balanced 
atomic honeycomb lattice resists the potential mismatch 
even after the cavity phase has reached its final value 
(p = 7r/4. It gets finally destabilized around r ~ 600 
and takes the form of a density-imbalanced honeycomb 
lattice, see Fig. [l6j As explained in Sec. |VB[ this struc¬ 
ture is not stable in the long-time limit and destabilizes 
subsequently into a triangular lattice as the atoms in the 
shallower sites move to the deeper sites. 
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FIG. 11: (Color online) Long-time dynamics of the system described in Fig. [t] (0 = 7t/2). Left panel: cavity photon numbers 
nj = \oi-j\ 2 (j = 1,2) as functions of time r = cjRt. As one can see, the cavity photon numbers fluctuate a lot and, in turn, 
so does the atomic honeycomb lattice. At the longest time r = 10 4 of our numerical simulation, the system has not yet fully 
reached its steady-state. Middle panel: corresponding contour plot of the atomic density. Right panel: effective lattice potential 
V e ff obtained at this longest time. The atomic honeycomb lattice is a stable structure but it takes much longer time to reach 
the steady-state. 



FIG. 12: (Color online) Transition from the density-imbalanced honeycomb atomic lattice to the triangular atomic lattice. 
Top-left panel: pump and cavity phase temporal sequences and time evolution of the cavity photon numbers. The cavity phase 
changes linearly from fa = ir/4: to (f>f=0 over the time interval At = 200. Meanwhile the pump strength U p decreases from 
0 to Uf — —20 over the time window 3Ar = 600. The other parameters are the same as in Fig. [4] Contour plots show the 
atomic density obtained at time r = ujRt — 200, before the cavity phase starts to decrease (top-right panel), at time r = 260, 
just before the triangular lattice become the dominant structure (bottom-left panel), and at time r = 400, when the cavity 
phase has reached its target value 4>f = 0 (bottom-right panel). The triangular lattice becomes the prominent structure around 
r ~ 280, well before the cavity phase reaches the target value 0. At this time, the cavity photon numbers rise considerably and 
then increase linearly in time like the pump strength for the rest of the time sequence. 
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FIG. 13: (Color online) Transition from the triangular atomic lattice to the density-balanced atomic honeycomb lattice. Top-left 
panel: pump and cavity phase temporal sequences and time evolution of the cavity photon numbers. The cavity phase changes 
linearly from fa — 0 (triangular lattice) to 0/ = 7r/2 (energy-balanced honeycomb lattice) over the time interval At = 200. 
Meanwhile the pump strength U p decreases from 0 to Uf = — 20 over the time window 3 At = 600. The other parameters 
are the same as in Fig. [4] When the cavity phase departs from 0, secondary minima start to grow and the effective potential 
becomes an energy-imbalanced honeycomb potential. However, the atomic density keeps its triangular lattice structure where 
atoms are sitting at the deeper minima of the effective potential, when the cavity phase reaches 7r/2 then the atomic density 
starts oscillating in time between the two triangular sub-lattices of the honeycomb lattice which now have the same depth. 
Top-right panel shows the atomic density at time t = 400 right at the time phase reaches ir/2. Bottom-left panel: atomic 
density distribution at time t = 402. The initial triangular lattice is still dominant. Bottom-right panel: atomic density 
distribution at time t = 404. The atomic distribution now lives mainly on the other triangular sub-lattice of the honeycomb 
lattice. 






















15 



0.8 

0.6 

0.4 

0.2 


- 2-1012 

xfk xfk 

FIG. 14: (Color online) Transition from the density-balanced atomic honeycomb lattice to the triangular atomic lattice. Top- 
left panel: pump and cavity phase temporal sequences and time evolution of the cavity photon numbers. The cavity phase 
changes linearly from 0i=7r/2to0/ = O over the time interval At = 200. Meanwhile the pump strength U p decreases 
from 0 to Uf — —20 over the time window 3Ar = 600. The other parameters are the same as in Fig. [4] Top-right panel: 
atomic density distribution at time r = 100, forming the expected density-balanced honeycomb lattice. As the cavity phase 
departs from 7t/2, consecutive potential wells start to have different depths. As long as the potential mismatch is weak enough, 
the density-balanced honeycomb lattice resists the phase change but, with increasing potential mismatch, finally gives in and 
becomes a fading density-imbalanced honeycomb lattice around r = 250 (bottom-left panel). When the cavity phase becomes 
even smaller, before it reaches 0 = 0 at r = 300, the atoms re-organize into the stable triangular structure. Bottom-right panel: 
atomic density distribution at time r = 600. 




VI. CONCLUSION 

In this work we have proposed to load a two- 
dimensional cloud of non-interacting cold bosonic atoms 
inside two identical initially-empty optical cavities with 
an angle of 27 t/ 3 between their axes. The atoms are 
driven by an external laser field. We have given the 
Hamiltonian of this hybrid system in the dispersive 
regime and we have derived the corresponding dynam¬ 
ical equations in the mean-field regime. The coherent 
superposition of the cavity and pump fields creates a dy¬ 
namical effective lattice potential with a triangular Bra- 
vais structure in which the atoms move. As a result of 
superradiance, the atoms self-organize into a triangular 
or a honeycomb lattice inside the cavities, the nature of 
the lattice depending on the relative phase between the 
pump and the cavity fields. Using the symmetry prop¬ 
erties of the effective potential, we have derived the dy¬ 
namical equations in reciprocal space and we have inves¬ 
tigated the condition required for superradiance to take 


place. Linear response theory shows that superradiance 
takes place irrespective of the relative phase between the 
cavity and pump fields. This is confirmed by our nu¬ 
merical results: atoms self-organize into a triangular lat¬ 
tice when pump and cavity fields are in phase, into a 
density-balanced honeycomb lattice when the fields are 
in quadrature and into a density-imbalanced honeycomb 
lattice in between. The stable atomic structures in the 
long-time limit are the triangular and density-balanced 
honeycomb lattices. The density-imbalanced honeycomb 
lattice only survives for a limited amount of time which 
becomes longer and longer as the cavity and pump fields 
approach quadrature. In the end, atoms redistribute 
equally into the triangular sub-lattice made of the deep¬ 
est potential wells. We have also studied the transition 
between these different lattice structures when the rel¬ 
ative phase between the fields is dynamically changed 
and confirmed the stability of the triangular and density- 
balanced honeycomb structures. A natural extension of 
this work would include the study of the survival time 
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FIG. 15: (Color online) Transition from the density-imbalanced to the density-balanced atomic honeycomb lattice. Top-left 
panel: pump and cavity phase temporal sequences and time evolution of the cavity photon numbers. The cavity phase increases 
linearly from fa = 7r/4 to <f)f = 7r/2 over the time interval At = 200. Meanwhile the pump strength U p decreases from 0 to 
Uf — —20 over the time window 3Ar = 600. The other parameters are the same as in Fig. Top-right panel: atomic density 
at time r = 200 where the cavity phase starts to increase. Later on, as the cavity phase increases, the shallower wells become 
deeper, however the atoms in deeper wells do not move to the shallower wells until phase becomes exactly 7r/2 and all sites 
have same depth. Bottom-left panel: atomic density at time r = 400 when the cavity phase reaches 7t/2 one gets an almost 
density-balanced atomic honeycomb lattice. From this time on, and for the rest of the time sequence, the atomic density starts 
oscillating fast in time between the two stable triangular sub-lattices. Bottom-right panel: hinting at this oscillating behavior, 
the atomic density at time r = 401 


of the density-imbalanced honeycomb lattices, the (pos¬ 
sibly chaotic) dynamics leading to superradiance when 
the fields are in quadrature and, importantly, the effect 
of atomic interactions on the self-organization process. 
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